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We calculate the expected galactic supernova neutrino signal at large next-generation underground 
detectors. At different epochs after the explosion, the primary fluxes can be quite different. For 
these primary neutrino fluxes, spectral splits induced by collective neutrino flavor transformations 
can arise for either mass hierarchy in both neutrino and antineutrino channels. We classify flux 
models according to the nature and number of these splits, and calculate the observable u e and 
Ve. spectra at Earth, taking into account subsequent matter effects. We find that some of the 
spectral splits could occur sufficiently close to the peak energies to produce significant distortions 
in the observable SN neutrino signal. The most striking signature of this effect would be presence 
of peculiar energy dependent modulations associated with Earth matter crossing, present only in 
portions of the SN neutrino energy spectra demarcated by spectral splits. These signatures at 
proposed large water Cherenkov, scintillation, and liquid Argon detectors could give hints about the 
primary SN neutrino fluxes, as well as on the neutrino mass hierarchy and the mixing angle #13. 

PACS numbers: 14.60. Pq, 97.60.Bw 



I. INTRODUCTION 

Observable effects of supernova (SN) neutrinos 1 in 
underground detectors are a subject of intense investi- 
gation in astroparticle physics. In particular, a lot of 
attention has been devoted to possible signatures of the 
Mikheycv-Smirnov-Wolfenstein (MSW) matter effect [l[ 
on neutrino flavor evolution [2tJ]- Shock- wave effects in 
the stellar envelope ((l-fl8j and the neutrino Earth matter 
crossing [l9l - |25| have been predicted to have dramatic 
effects on SN neutrino oscillations and to produce signa- 
tures that could allow us to extract valuable information 
on the unknown neutrino mass and mixing parameters, 
like the neutrino mass hierarchy and the mixing angle #13. 

In the recent years, it has been understood that the 
paradigm of neutrino flavor transformations in super- 
novae, based on the only MSW effect with the ordinary 
matter [l9l ]. was incomplete. Novel phenomena have 
been found to be important in the region close to the 
ncutrinosphcrc, where the neutrino density is so high 
that effects of neutrino-neutrino interactions dominate 
the flavor evolution. These effects have been understood 
and characterized in an increasingly realistic way in a 
long series of papers [2fj| - t73j . 

The neutrino-neutrino interactions provide a large po- 
tential for the neutrinos due to the neutrinos themselves, 
which causes large and rapid conversions between dif- 
ferent flavors. The transitions occur collectively, i.e. in 
a coherent fashion, over the entire energy range. Unlike 



1 Neutrinos will refer to both neutrinos iy) and antineutrinos (y) 
in this article. 



ordinary neutrino oscillations, these collective oscillations 
depend strongly on the original SN neutrino fluxes. In a 
supernova v e and v e are distinguished from other flavors 
due to their additional charged-current interactions. The 
v^ T and their antiparticles, on the other hand, are 
produced at practically identical rates. Following the 
standard terminology, we define the two relevant non- 
electron flavor states as v x . y = cos 023^T sm #23^V , where 
623 — 7r/4 is the atmospheric mixing angle. Since the 
initial v x and v y fluxes are identical, the primary neutrino 
fluxes are best expressed in terms of v e , v e and v x . The 
collective pair conversions v e v e O v xy v xy take place 
even for extremely small mixing angles within the first 
O(500) km 34 1, much before the MSW flavor conversions 
start. In a typical supernova, enough asymmetry between 
neutrino and antineutrino number fluxes is expected, 
thus complicated angle-dependent decoherence effects 
due to the anisotropic neutrino emission in the SN 
environment are likely to be suppressed [13] ■ There may 
be however a slow-down of collective flavor conversions 
due to multi-angle effects in deep supernova regions [72| . 

During the accretion phase (t < 0.5 ms after the 
core-bounce) in typical supernova simulations [74T - t78j 
one finds an almost perfectly equipartitioned luminosity 
between the neutrino flavors: L„ e ps L Pe « L Ux . This has 
been the benchmark scenario in most studies of collective 
flavor oscillations, which predicts: For normal mass 
hierarchy (NH) v and T> remain unaffected by collective 
oscillations. For Inverted hierarchy (IH), the end of 
collective oscillations is marked by an almost complete 
exchange of the v e and v y flavor spectra for u's, while for 
v's the swap occurs only above a characteristic energy 
fixed by lepton number conservation, giving rise to a 
spectral split in the v energy distributions [H, SO, H3 • 
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It is not obvious that the neutrino fluxes maintain 
their approximate cquipartition even at late times. In- 
deed, the dclcptonization of the core is probably faster 
than the cooling of the proto-neutron star, so that the 
asymmetry between the fluxes can become smaller and 
the spectra at late times depend on details of the neutrino 
transport (the density and the temperature profiles, as 
well as the treatment of the interaction rates). Due 
to this complexity, results from SN simulations are not 
unambiguous. Late-time cooling calculations performed 
by the Garching group [76| show a cross-over of the 
different neutrino luminosities L Ua , > L Ve L Pc , i.e. 
at late times the v x flux becomes relatively larger [77j . 
Similar results, but less pronounced, have also been 
obtained recently by the Basel group 78 | and in the most 
recent long-time simulations of the Garching group [79j . 

This latter case (L Vx > L Vc w i Pc ) was recently 
studied [HI, [7(1 (zH , and one found the occurrence of un- 
expected multiple spectral swaps and consequent spectral 
splits for both j/'s and v's in either mass hierarchy. The 
lesson from this result is that the benchmark scenario 
(one spectral split in v e spectrum and complete swap 
in the D e for inverted hierarchy) is in fact a special 
case, while the phenomenology of the spectral features 
can be more complex. To this end, a detailed study 
has been performed in [63}, where the impact of the 
variations of the neutrino luminosities has been explored, 
finding abrupt changes in the number and the position 
of the splits by slightly changing the ratio between the 
luminosities of the different neutrino species around some 
critical value. 

Given the sensitivity of the self-induced spectral splits 
on the original SN neutrino emission features, their 
detection in a galactic SN neutrino burst could provide 
a tool to reconstruct the original neutrino fluxes. Since 
spectral splits could potentially affect both neutrinos and 
antineutrinos, it seems worthwhile to investigate comple- 
mentary detection techniques with sensitivity to v e and 
v e respectively. Around the world, there is an active 
discussion about the feasibility of three different classes 
of large underground detectors for low-energy neutrino 
physics and astrophysics, viz. water Chcrcnkov detectors 
with fiducial masses ranging from 0(0.4 Mt) [80l - [8^ | to 
5 Mt [HI, liquid scintillation detectors with masses of 
0(50 kt) [H] and liquid Argon Time Projection Cham- 
bers (LAr TPC) with fiducial masses of 0(100 kt) [H- 
HH. In particular, these three detection techniques are 
the backbones of the European project LAGUNA (Large 
Apparati for Grand Unification and Neutrino Astro- 
physics) [8i| and the LBNE (Long Baseline Neutrino 
Experiment) in DUSEL (Deep Underground Science and 
Technology Laboratory) [90] . The physics potential 
of such devices for supernova neutrino detection would 
be extremely high. In particular, water Cherenkov and 
scintillation neutrino experiments are mostly sensitive 
to supernova v through the inverse beta decay process 
D e +P —> e + +n. On the other hand, LAr TPC would have 
a high sensitivity to SN v, through the charged current 



interactions of v e with the Ar nuclei in the detector. 
The complementarity between these different detection 
techniques would allow to compare the different features, 
like the spectral splits, occurring in SN neutrino spectra. 

Motivated by this intriguing perspective, our aim in 
this paper is to calculate the observable SN signal at 
these detectors. The main feature of our work is an im- 
proved treatment of the oscillation physics, particularly 
collective effects, in different phases of the SN explosion. 
We also investigate the detectability of possible features 
in the observable galactic supernova signal at future 
large underground detectors and examine the physics 
potential. 

The plan of the article is as follows. In Sec. [n] 
we describe the parameter space of expected primary 
neutrino fluxes, the equations for the neutrino flavor evo- 
lution, and neutrino mixing. In Sec. IIII1 we give a short 
overview of the collective effects on neutrino spectra and 
the multiple spectral splits. We also perform a scan on 
the SN flux parameter space to identify different classes 
of fluxes that lead to qualitatively different final spectra. 
Sec. IIVI discusses the subsequent flavor conversions due 
to MSW effects and possible Earth matter effects, and 
presents the net survival probabilities for v e and v & 
arriving at the detectors. In Sec. [V] we describe the 
features of our reference detectors (fiducial mass, cross 
sections and energy resolution). In Sec. |VT] we present 
our results on the detectability of the spectral signatures 
in the v e and v e signal. In particular, we discuss the 
distinctive features due to interplay of Earth matter 
effects and spectral splits. Wc also briefly comment on 
neutrino oscillation effects on the early ncutronization 
burst and in relation to the SN shock-wave propagation. 
Finally, in Sec. lVIll we comment about our results, discuss 
future perspectives, and conclude. 



II. FLUXES, POTENTIALS AND MIXING 
PARAMETERS 



In this section, we describe the primary fluxes, the 
effective potentials they encounter during their propa- 
gation, and the neutrino mixing parameters that we use 
in our numerical simulations. 



A. Primary neutrino fluxes 

The primary neutrino fluxes of all neutrino and an- 
tincutrino species can be written as 



(1) 



where v = v e ,V e ,v x , ip(E) is the normalized primary 
neutrino spectrum, i.e. J dE ip(E) = 1, and is 
the total number flux. For <p(E) we use the spectral 
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parameterization given in Rcf. [76j |: 

K + ( E \ a » e 



tp(E) 



T(a u + 1) 



(2) 

where (E v ) is the average v energy, ol v is a spectral pa- 
rameter, and r is the Euler gamma function. The values 
of the parameters are model dependent [74| [75| . There is 
significant variation across different SN simulations, but 
we find that the (E Vt ,) and (E Pb ) are usually close to 

(E„ e ) = 12 MeV and {E Pe ) = 15 MeV , (3) 

whereas {E V:c ) tends to vary in the range 

16 MeV < (E U J < 25 MeV . (4) 

We fix the spectral pinching parameter at a„ = 3 for all 
species. 

The total number fluxes $° are related to the lumi- 
nosities of the neutrino species through 



L, 



{E v ) 



(5) 



For definitcness, we assume a typical supernova which 
emits Eb « 3 x 10 53 erg in a duration of t = 10 s. The 
time integrated luminosities of all species combined is 
equal to the total emitted energy: 



dt(L Ue +L Pe + 4L„J =E E 



(0) 



such that the average number fluxes are about 
10 55 s _1 MeV _1 , and about 10 58 neutrinos are emitted 
in the entire SN explosion. 

The ratio between the different luminosities varies 
significantly across models, but we find that (See Fig. 3 
of [13 and Table 7.3 of (Hi for example) 



E Vc I L Ve ~ 1 , 
0.5 < L v jL Ve < 2.0 



(7) 



The relative number fluxes change over the duration 
of the neutrino burst. In particular, the fluxes in the 
accretion phase are not necessarily similar to those in the 
cooling phase. These primary fluxes are subject to flavor 
conversions, of which the collective oscillations that we 
will discuss in the next section depend quite sensitively 
on the fluxes themselves. 



B. Neutrino equations of motion 

Mixed neutrinos are described by matrices of density 
p p and p p for each (anti)ncutrino mode. The diagonal 
entries are the usual occupation numbers whereas the off- 
diagonal terms encode phase information. The equations 
of motion (EoMs) are [92| 



idtpp = [Hp,p p ] 



where the Hamiltonian is 



Hp — ftp + Vmsw + V„„ 



(9) 



Here fi p is the matrix of the vacuum oscillation frequen- 
cies for neutrinos, which is 



Q p = diag(m 1 ,m 2 ,m 3 )/2|p 



(10) 



in the mass basis. For antineutrinos fi p — > — fi p . The 
matter effect due to the background electron density n e 
is, in the weak interaction basis, i.e. (v ei v^, v T ), 



Vmsw = V / 2G i ?n e diag(l, 0, 0) , 



(11) 



neglecting the second-order difference between the 
and v T refractive indices that could be important for 
collective neutrino oscillations only at really early post- 
bounce times it < 300 ms) Hi| [92| . 

The effective potential due to the neutrino-neutrino 
interactions is 



V„ V = V2G F J ^S_( Pq -p q )(l 



,) , (12) 



where v p is the velocity of the neutrino mode p. 

In spherical symmetry the EoMs can be expressed as 
a closed set of differential equations along the radial 
direction (3?i [H|. The factor (1 — v q ■ v q ) in the 
Hamiltonian gives rise to "multi-angle" effects for neu- 
trinos moving on different trajectories [Hj]. However, 
for realistic supernova conditions the modifications are 
presumably small (except for a shift in the onset of 
the flavor conversions [72[)j allowing for a single- angle 
approximation (37l . [60j . We implement this approxima- 
tion by assuming all neutrinos to have been launched 
from the neutrinosphcrc at radius R, with an angle 45° 
relative to the radial direction [37], In this case, 

the formal derivation of the single-angle approximation 
explicitly gives the radial dependence of the neutrino- 
neutrino interaction strength as [53| 



R 

y vv {r) oc p{r) = po — Q 
Here, po is the interaction strength 



Mo 



= V2G F (F°-Fl) , 



(13) 



(14) 



the r 2 scaling comes from the geometrical flux dilution, 
while the collinearity factor 



1 - y/1 - {R/rf 



2 



(R/r) 



■Y2 



1 /R 

2 W 



for r — > oo 



(15) 



arises from the (1 — v p ■ v q ) structure of the neutrino- 
neutrino interaction. The decline of the neutrino- 



(8) neutrino interaction strength, p(r) 



for r > R, 
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is evident, while the numerical coefficient for large r 
depends on the launching angle, here taken to be 45°. 
When fluxes at the ncutrinosphere radius R = 10 km are 
taken, we get 



/i = 7 x 10 5 km" 1 



(16) 



Of course, the physical neutrinosphere is not a well- 
defined concept. Therefore, the radius R simply rep- 
resents the location where we fix the inner boundary 
conditions. However, essentially nothing happens close to 
the neutrinosphere because the in-medium mixing angle 
is extremely small. Therefore, as far as the vacuum 
and matter oscillation terms are concerned, it is almost 
irrelevant where we fix the inner boundary condition. 



C. Neutrino mixing parameters 



We work in the rotated basis |47| 

(v e ,v x ,v y ) = R T {923,){v e ,v ll ,v T ) . 

This is equivalent to taking #23 = in the neutrino 
mixing matrix, which makes no difference to v e and v e 
evolution as long as the primary fluxes of and v T are 
identical, as we have assumed. The values of the other 
mixing angles are sin 2 B\2 — 0-30 and sin 2 #13 < 0.04 (§3] 
in vacuum. The neutrino mass-squared differences in 
vacuum are taken to be Am 2 tm = 2 x 10 -3 cV 2 and 
Am 2 ol = 7.6 x 10 -5 eV 2 , close to their current best-fit 
values [94| . We study both the cases of normal neutrino 
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'atm 
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> 0) and 



mass hierarchy (NH: Am 
inverted mass hierarchy (IH: Am 
We ignore possible subleading CP violating effects [5f 
by setting <5cp = 0. 

Matter effects in the region of collective oscillations (up 
to a few 100 km) suppress the mixing angles and slightly 
modify the neutrino mass-square differences. We take 
the matter-suppressed mixing angles to be #13 
10~ 3 and the effective mass-square differences Am 2 



^12 
"^atm 



Ar 
0.4A 



COS 01 



and At 



'sol 



Am 2 j cos 6<i2 



'"I 



sol 



|54(. Apart from these shifts, matter 
effects typically do not disturb the development of the 
collective neutrino oscillations, except at very early times 
it < 300 ms) when the effective electron density n e would 
become larger than the neutrino density n„, suppressing 
the collective flavor conversions [54 ] ■ For simplicity, in 
the following we will always focus on later times where 
collective oscillations are not inhibited by a strong matter 
term. MSW conversions occur after collective effects have 
ceased [13, [Ml ■ Therefore, their effects factorize, and 
will be described in Sec. IIVI 



III. COLLECTIVE FLAVOR CONVERSIONS 

Although the subject of collective oscillations has been 
discussed in great detail in number of papers, most have 



focused on the benchmark scenario used in [42j. It was 
pointed out in [6^, [zO, [zj that the final outcome of 
collective effects can be quite complicated, depending 
on the primary fluxes. In this section we will perform 
a detailed study investigating the effect of the primary 
fluxes on flavor oscillations during the different post- 
bounce phases. 



A. Spectral swaps and splits 

Near the neutrinosphere, due to the large neutrino 
density, the neutrino-neutrino interaction energy is very 
large. This ensures that neutrinos exhibit collective 
oscillations, i.e. neutrinos of all energies oscillate coher- 
ently with the same average frequency. As the neutri- 
nos stream outward and the neutrino density becomes 
smaller, self-induced oscillations begin to take place 34 1. 
These oscillations, being instability driven, have a large 
amplitude even for a very small mixing angle. In the 
presence of a slowly decreasing background neutrino 
density, the collective effects on the neutrinos become 
negligible beyond a point, and their imprint on the 
spectra is left in the form of the spectral splits. The 
energy and the number of the splits is crucially dependent 
on the neutrino mass hierarchy (normal or inverted), and 
on the relative sizes of the primary number fluxes of 
different flavors. 

The resultant dynamics of the spectral splits has been 
elaborated in [62l [7l| : spectral swaps can develop around 
unstable crossing points. Here crossing points are those 
energies E c where spectra of different flavors cross each 
other, i.e. 



F? e (E c ) = Fl(E c ) or F° (E c ) = i£ (E c ) . (17) 

A given crossing point is unstable for e y ("atmo- 
spheric sector") swap, triggered by (Am 2 tm ,#i3), if at 
the energy E c , 



d(F° e - F° x )/dE < for IH , 
d(F? e -F?J/dE>0 forNH 



(18) 



and analogously for antineutrinos. Thus, in IH we get 
a v e f-> Vy or v e v y swap around every crossing 
with a negative slope, while in NH we get such a swap 
around every crossing with a positive slope. Each swap 
is demarcated by two spectral splits. Two subtle points 
may be noted: (i) If two unstable crossings are very close 
to each other, the swaps may be influenced by each other 
and even merge, and (ii) If an unstable crossing is too 
closely flanked by stable crossings, the width of the swap 
around it is suppressed exponentially, which also leads 
to adiabaticity violation leading to smearing out of that 
swap [62j . 

In addition, a crossing point is unstable for the e x 
("solar sector") swap, triggered by (Am 2 ol , 6*i 2 ), if [7l| 



d(F° -F°)/dE>0 



(19) 
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FIG. 1: Regions of the flux parameter space that produce qualitatively different spectra after collective flavor transformations. 
Left panel for NH, Right for IH. The bands represent the positions of spectral swaps in energy, and their nature: vertical stripes 
correspond to u e o v v swap, horizontal stripes correspond to v e f-> v x swap, unstriped regions have no swap. The widths of the 
striped / unstriped regions approximately indicate the energy range where they are present, the complete band corresponding 
to the energy range 1-50 MeV. 



and analogously for antineutrinos. We thus obtain a 
v P -H- v x or v e f-> v x swap around every crossing with 
a positive slope. These swaps tend to be more non- 
adiabatic than those in the atmospheric sector, since 
the natural frequency of the solar sector, Am^ ol /(2E), is 
smaller than that of the atmospheric sector. As a result, 
the rate of change of the collective interaction strength 
may become too fast for there to be enough time for the 
solar sector instability to grow. 

The combination of the atmospheric and solar sector 
swaps may give rise to three-flavor effects, pointed out 
recently in [7l| . For NH, mass splittings are positive 
in both the atmospheric and solar sectors, and both 
the mixing angles are small due to matter suppression. 
Both the sectors then tend to produce instabilities at the 
same energies, but the atmospheric sector wins due to its 
higher natural frequency (and hence larger adiabaticity). 
As a result, only v & -H> v y and v & <-t v y swaps develop, 
and no three-flavor effects are observed. On the other 
hand, in IH the atmospheric and solar instabilities are 
in different parts of the spectrum. The solar sector in- 
stability then operates in the high-energy region without 
hindrance and causes an additional swap that may even 
merge with the atmospheric sector swap, partially erasing 
one of the spectral splits. However, the split usually is 
not completely erased because of the non-adiabaticity of 
the swap in the solar sector. 



B. Flux dependence of spectral splits 

The number and positions of the spectral splits de- 
pends on the primary neutrino spectra [63| . However a 
clear mapping between primary fluxes and split patterns 
is still lacking. The flux parameters with the widest 



variation among models are the average energy of the 
v x flux (-Ei, x ), and the v x luminosity L Vx . We scan the 
SN neutrino parameter space in the variables L Uas /L Vc 
and {E V:c ) in the ranges defined by Eqs. (0]) and (JT]), and 
numerically evolve the system of neutrinos to determine 
the regions of the flux parameter space where different 
spectral swaps are effective. We then classify the primary 
spectra depending on the number and nature of these 
swaps. 

The results of this classification may be represented 
by a "phase diagram" shown in Fig. [T] Note that the 
idea of a phase diagram, in which the oscillation patterns 
change a bruptly across regions, has been implemented 
earlier in [63, [671 where a two-flavor analysis of collective 
effects is carried out. Our survey involves a three-flavor 
treatment of collective effects, and is carried out over a 
different slice of the parameter space, i.e. we allow the 
v x ener gies to vary while fixing Lp e = L Ve . The surveys 
in Ref. [63l . [67| chose to fix the v x energy and vary the 
relative fluxes / luminosities of the flavors. Our choice is 
motivated by the observation that most simulations tend 
to predict L Ve w i Pc , whereas (E Vx ) turns out to be much 
less robustly predicted. We therefore believe that this 
survey provides useful information that is complementary 
to the one in [63|,|67|. 

From the phase diagram, it is observed that for NH, 
there are two phases with qualitatively different split 
patterns. Phase A, which approximately corresponds to 
L Ve ^ L Ux , shows no spectral swaps in v or v. On the 
other hand, in the other phase C with L Ve < L Vx , one 
obtains a v e <-> v y swap and a v e <-> v y swap, both at 
energies > 25 MeV. The situation is more complicated 
for IH. Here the phase A, which is approximately the 
same as the phase A in NH, corresponds to a v e O v y 
swap at energies > 8 MeV, and a v e <H- v y swap at 
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FIG. 2: Fluxes (integrated over 4-7T solid angle) for the benchmark models: A (left panel), Ci (right panel). In each panel, 
fluxes are shown for antineutrinos (left) and neutrinos (right) at initial time at neutrinosphere (top), post-collective oscillations 
in IH (middle) and NH (bottom). The e flavor is shown in red/dotted, x in green/solid, and y in blue/dot-dashed. The fluxes 
for models C2, C3, C4 are not shown, as they are similar to Ci, with certain features degraded (see text). 



energies > 4 MeV. The phase C however has different 
features depending on its sub-phases Ci, C2, C3 and C4. 

• Ci correponds to a large L Vsc and large (E V J), and 
hence here the spectral swaps are more pronounced. 
Both v e <H- v y as well as v e v y swaps appear at 
intermediate energies 10 MeV < E < 25 MeV. For 
E > 25 MeV, there are v e -H> v x and v e o v x swaps. 

• In C2, the energy {E U J) is lower as compared to 
Ci. Then v e and v y spectra are very similar at 
intermediate energies. As a result, the v e -f-> v y 
swap at intermediate energies has low adiabaticity, 
and is not effective. 

• In C3, the luminosity L„ x is lower as compared to 
Ci. Then the v e and v x spectra, as well as the v e 
and v x spectra, are very similar at high energies. As 
a result, the v e <-> v x and v e f-> v x swaps at high 
energies have low adiabaticity, and are not effective. 

• C4 corresponds to almost equal v e and v x number 
fluxes everywhere except at intermediate energies 
near the peaks. As a result, the v e <-> v x and v e <-> 
v x swaps at high energies are non-adiabatic and 
hence ineffective. The v e <-> v y and T> e <-> T> y swaps 
at intermediate energies are partially adiabatic. 
This case can be thought to be the combination 
of C2 and C3, with only the common features 
surviving. 



In Fig. [3 we illustrate the spectral swaps for two 
benchmark points 

A : (E Vm ) = 18 MeV, L V JL V , = 0.7 , 
Ci : (E Vm ) = 20 MeV, L V JL V , = 1.5 . (20) 

The benchmark point A lies close to the predictions of the 
Garching flux model [75| during accretion phase, while Ci 
lies close to the predictions during cooling phase. Here 
we show the primary fluxes of all v as well as v species, 
and show how they change after collective effects. The 
v e and 77 e spectra after collective oscillations exhibit the 
spectral features described before. We point out that 
for the Ci flux, in IH the oscillated v e spectrum shows 
a single split, the one at low-energy, corresponding to 
the swap with v y , while the high-energy split is canceled 
by the further v e f-> v x swap. However, this high- 
energy split is still visible in the v x and v y spectra. 
We will see in the following that this high-energy split 
would reappear in the final v e spectrum at Earth, since 
that is a superposition of the different v spectra after 
collective oscillations. The same observations apply for 
the v e spectrum. Moving from the case Ci to C4 in 
the phase C one may observe a gradual degradation in 
the observability of the different spectral splits in final 
spectra (not shown). 
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IV. FLUXES AFTER MSW CONVERSIONS 
A. Survival probabilities 

After collective oscillation die out, the primary neu- 
trino fluxes F® experience spectral swaps and then 
further undergo the traditional MSW conversions in 
SN [9{| while passing through the resonance regions 
H and L, corresponding to flavor transitions in the 
two-neutrino sectors associated with (Am 2 tm ,0i3) an d 
(Am 2 ol , #12), respectively. After exiting the resonance 
regions, the neutrino mass eigenstates travel indepen- 
dently until they reach Earth, wherein they are detected 
as flavor eigenstates. The fluxes of v e and v e arriving at 
Earth can be written as uM 



Fu c = pF° e + (l-p)F£ 
Fv e = PFI + (1-P)F? 



(21) 



where p and p are the v e and v e survival probabilities, 
respectively. The survival probabilities are determined 
by the adiabaticity of the MSW £f-resonance and in 
general are sensitive to the neutrino energy and to the 
SN matter density profile (see, e.g., Q). This dependence 
vanishes for a large time-window in the limiting cases of 
adiabatic transitions (sin 2 6*13 > 10 -3 ) and strongly non- 
adiabatic transitions (sin 2 #13 < 10~ 5 ) (see, e.g., [Hj])- In 
the following, we shall consider for simplicity only these 
two limiting cases. 

In these two regions of the #13 parameter space, the 
expressions for p and p can be written down in a simple 
form, as shown in Tables U and fll] The survival prob- 
abilities have been given separately for the low-energy 
(E < E\ ow ), intermediate energy (£q ow < E < -Ehigh) 
and high-energy (E > i?high) regions. Here E\ ow and 
-Ehigh, respectively, are the energies where the low-energy 
and high-energy spectral splits are likely to be present. 
For v ei one has £i ow w 10 McV and E high ~ 25 MeV. 
For v e , one gets E\ ow ps 5 MeV and E high w 25 MeV. 
The survival probabilities also depend on the primary 
fluxes, and are given for the phases A (L Ve > L Vm ) and 
C {L Ue < L Vx ), covering the complete parameter space. 
Some special subcases of C are given in the caption. 
These are related to the cases C2, C3, C4 shown in Fig. 1, 
for which we have seen that some of the splits can be 
suppressed by the presence of too close v e and v x spectra 
or by the low v x luminosity. 

Note that the survival probability is never 1. The 
vacuum mixing due to the angle Q12 always mixes the 
v e or v e spectra with the other flavors. Effects of the 
spectral swaps are invariably reduced by a factor of 
sin 2 12 ~ 0.3 or cos 2 12 « 0.7. 



B. Earth matter effect 

The neutrino survival probabilities at the Earth given 
in Tables HI and [TTI are calculated assuming that neutrinos 



escaping the star travel through vacuum before reaching 
the detector. If the supernova is shadowed by the 
Earth for a detector, the neutrinos will travel a certain 
distance through the Earth, and will undergo Earth 
matter oscillations during this propagation. Since the 
neutrinos arrive at the Earth as mass eigenstates, the net 
effect of these oscillations can be written in terms of the 
conversion probabilities Pi e = P{i>i — > v e ). Neglecting 
the effect due to small #13, one can obtain the net survival 
probabilities by the substitution [l9| 

(cos 2 flu, sin 2 12 ) -> (1 - P 2e , P2e) (22) 

for neutrinos in Table U and 

(COS 2 012, sin 2 012) -> (1 - P 2e , Pie) (23) 

for antineutrinos in Table HB Here 

P 2e = P{v 2 ^v e ) , and P 2e = P{D 2 ^D e ) (24) 

while propagating through the Earth. Analytical ex- 
pressions for P 2e and P 2e can be calculated for the 
approximate two-density model of the Earth [23| . When 
neutrinos traverse a distance L through only the man- 
tle of the Earth, these quantities have a very simple 
form El, [II: 



P 2e = sin 2 0i2 + sin 20^ x 

sin(20»-20 12 )sin 2 ( L 



AE sin 20; 



P 2e = sin 2 0i2 + sin 20^ x 



(25) 



(26) 



sin(20™-20 12 )sin 2 '' Am -' Sin ^ 12 



4E sin 20 



L 



12 



where 0^ and 0^ are the effective values of 0i2 in Earth 
matter for neutrinos and antineutrinos, respectively [95| . 

The Earth crossing thus induces an oscillatory signa- 
ture in the neutrino energy spectrum. However note 
that these oscillations are absent if the corresponding 
survival probabilities vanish: they occur only when p =/= 
or p ^ in Tables Q] and [TTJ Depending on the swap 
pattern, the Earth effect can then appear in different 
energy regions in the spectra. As we shall see in Sec. I VII 
the mere observation of Earth matter effects is sometimes 
enough to distinguish between different flux and mixing 
scenarios. 



V. NEUTRINO DETECTION 

In this section we describe the main aspects and 
ingredients of our calculations of supernova neutrino 
event rates. The oscillated SN neutrino fluxes at Earth, 
F v , must be convolved with the differential cross section 
<7 e for electron or positron production, as well as with 
the energy resolution function R e of the detector, and 
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TABLE I: Survival probability p for u e at low, intermediate and high energies, for fluxes in phase A and C. Within phase C, 
the exceptions to the rule are denoted by the brackets (..) for C3 and C4. 
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TABLE II: Survival probability p for v e at low, intermediate and high energies, for fluxes in phases A and C. Within phase C, 
the exceptions to the rule are denoted by: brackets (..) for C3 and C4, square brackets [..] for C2 and C4. 





Phase A (L„ e > L„J 
E < Blow Bi „ < B < Bhigh B > Bhigh 


Phase C {L Ve > L v J 
B < Bio„ Blow < B < Bhigh B > Bhigh 


NH sin 2 0i3 > 10" 3 
sin 2 0i3 < 10" 5 


COS 2 012 COS 2 012 COS 2 012 
COS 2 012 COS 2 012 COS 2 012 


COS 2 012 COS 2 012 
COS 2 012 COS 2 012 


TTJ sin 2 0i3 > KT 3 
sin 2 0i3 < 10" 8 


COS 2 012 COS 2 012 
COS 2 012 


cos 2 0i2 [0] sin 2 0i2 (0) 

COS 2 012 [COS 2 012] Sin 2 012 (cos 2 012) 



the efficiency e, in order to finally get observable event 



rates 
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N, = F v 



(27) 



To calculate the total number of events, we will assume 
the supernova distance d = 10 kpc and integrate the 
event rates over t = 10 s, assuming the fluxes to 
be constant over the entire duration. The fluxes are 
thus to be thought of as the time-averaged fluxes over 
10 s. Alternately a more detailed modeling of the time- 
dependent flux is required, in which case a time-binned 
analysis may be performed. 

We will now describe the main characteristics of three 
types of detectors: water Cherenkov detectors, scin- 
tillation detectors, and liquid Argon Time Projection 
Chambers, that we have used to calculate the signals. 



A. Water Cherenkov detectors 

In large water- Cherenkov detectors, the golden channel 
for supernova neutrino detection is the inverse beta decay 
of electron antincutrinos 2 



P 



(28) 



For this process, we take the differential cross section 
from [9t| . The total cross section grows approximatively 



2 We will neglect the sublcading neutrino interaction channels in 
the detectors, assuming that they can be separated at least on a 
statistical basis. 



as E 2 . Wc fold the differential cross sections for e + 
production with a Gaussian energy resolution function 
of width A. The value of A is predominantly determined 
by the photocathode coverage of the detector. For our 
calculations we assume 0, 



A wc /MeV = 0.47V^e/MeV , 



(29) 



where E e is the true positron energy. A galactic SN is 
expected to produce 0(1O 5 ) events in a Mt-class water 
Cherenkov detector with a fiducial volume of about 
400 kt [H. 



B. Scintillation detectors 

In liquid scintillators, the main channel for SN neutrino 
detection is the inverse beta decay of £Vs, the same as 
that in water Cherenkov. However, here the positrons are 
detected through photons produced in the scintillation 
material. Since a larger number of photons can be 
produced in a scintillation detector, these have typically a 
much better energy resolution than the water Cherenkov 
detectors. Indeed, the energy resolution of a scintillation 
detector may be better by more than a factor of 6. 
The energy resolution of the scintillator detectors is 
determined by the number of photo-electrons produced 
per MeV, which for this type of detectors is expected to 
be given by as good as [9£| [9fj| 



A S c/MeV = 0.075 v'^eyMeV 



(30) 



Since the Earth matter oscillations described in the 
previous section may get smeared out by the finite energy 
resolution of the detector, it is clear that the energy 
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resolution plays a crucial role in the efficiency of detecting 
Earth effects. For a fiducial mass of 50 kt, one expects 
O(10 4 ) events for a galactic SN || . 



C. Liquid Argon Time Projection Chambers 

LAr TPC detectors would be particularly sensitive 
to SN electron neutrinos through their charged current 
interactions with Argon nuclei 

u e + 40 Ar -> 40 A'* + e~ , (31) 

which proceed via the creation of an excited state of 
40 A and its subsequent gamma decay. The Q-value for 
this inverse beta decay process is 1.505 MeV. The cross- 
section for this charged current reaction is taken from 
[lOOj ] . Energy resolution in such detectors is expected to 
be very good. The energy resolution for leptons in liquid 
argon time projection chambers has b een c alculated by 
the ICARUS collaboration who report jl0l| 

A LAr /MeV = 0.11 v/Ee/MeV + 0.02 E e /MeV . (32) 

The fiducial volu me fo r supernova neutrino detection is 
taken to be 100 kt [l02j |. With these standard inputs, one 
expects C(10 4 ) events from the interactions considered 
above (89|, while the event rates produced by the other 
interaction channels are smaller by at least an order of 
magnitude. Due to the strong sensitivity to v e , the liquid 
Argon technique would be complementary to the water 
Cherenkov and scintillation detectors, which are mostly 
sensitive to i> e 's. 



VI. OBSERVABLES SENSITIVE TO NEUTRINO 
FLAVOR CONVERSIONS 

In this Section we discuss about possible signatures 
of SN neutrinos flavor conversions observable in the 
neutrino energy spectra in the three different types 
of detectors, presented before. First we calculate the 
neutrino energy spectra at different detectors and we 
discuss about the observability of the spectral splits and 
of the Earth matter effect. Then, we also comment about 
two other possible observables sensitive to the neutrino 
mixing, namely the v e prompt neutronization burst and 
the shock- wave effect on the SN neutrino signal. 

In order to get a realistic idea of what features of the 
neutrino spectra can be observable at a neutrino detector, 
we show the events observable at different detectors for 
two of our benchmark cases, A and Ci, in Figures [3] 
and EJ respectively. For the v e spectra, we show 
the events at a LAr TPC , while for the v e spectra, we 
use the scintillation detector. The events spectra at the 
water Cherenkov can be obtained by smearing the v e 
spectra at the scintillation detector. We do not show 
them separately in this figure. The event-rates will scale 



as Eb/cP and the volume of the detector, and one can 
read off the statistical errors, that go as y/N, accordingly. 

We now analyze features in the detected neutrino spec- 
tra that can tell us about the flux and mixing scenarios, 
and examine the feasibility of detection of these features 
under different conditions. 



A. Spectral splits 

Spectral splits are produced at the boundaries of 
swaps, as shown in Fig. [21 At a spectral split, the survival 
probability p or p jumps suddenly, so that a sharp jump 
in the spectrum is in principle observable if the primary 
spectra of electron and non-electron flavors at the split 
energy are different. The magnitude of the jump is 
reduced by the nonzero value of #12, and its sharpness is 
somewhat smeared out by the energy resolution as well 
as possible non-adiabatic nature of the swap. 

Figures [3] and |4] bring out the following features of the 
spectral split: 

• In the phase A, collective oscillations occur in IH, 
producing a spectral split in the z/s and a complete 
swap of i>'s. However, the split energy is relatively 
low (around 10 MeV). At this energy, the difference 
in F v ^ and F v=c is small, and the charged current 
Verhr cross section, being cx E 2 , is also small. 
As a result, this "classic" spectral split in phase 
A is practically unobservable, as can be seen in 
Fig. [3] We realize that the peculiar interplay 
between collective and MSW effects can lead to 
quite different neutrino spectra for the two mass 
hierarchies and for large and small values of 6*13. 
However, lacking of a calibration of the original SN 
neutrino fluxes, in this case it would be hard to 
make any strong statement about neutrino mixing 
just from the observation of the neutrino energy 
spectra. 

• Observable spectra produced by v with primary 
fluxes as in Ci are represented in Fig. 01 In the 
phase C, splits are possible in both, v e and v e 
spectra at higher energies, around E sp \a = 25 MeV. 
This is near the peak of the primary spectra, 
so the difference between the electron and non- 
electron flavor spectra is more pronounced. This 
makes possible the detection of signatures related 
to spectral splits. In particular for 77 e , in NH and 
in IH with sin 2 6*13 > 10~ 3 , the observable positron 
spectrum in a liquid scintillation detector (or in a 
water Cherenkov detector) is mostly due to F£ for 
E < E sp \it and at higher energies (see Table II) . 
This would produce a bimodal positron spectrum, 
with two peaks corresponding to the two peaks 
of the initial antincutrino distributions. Instead, 
at sin 2 #13 < 10~ 5 in IH, the positron spectrum 
above E ~ 10 MeV will be mostly produced by F£ 
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FIG. 3: For benchmark flux model A, v e and v R energy spectra at 50 kt scintillator and 100 kt LAr TPC detectors, in both 
the hierarchies NH (upper panels) and IH (lower panels), with and without oscillations due to Earth matter effects. The 
spectra with Earth matter effects (EM) have been calculated for L = 8000 km through the Earth, and have been denoted 
by thinner lines. The spectra at a water Cherenkov detector can be obtained by smearing the energy of v e at the scintillator 
detector, and multiply the number of events ~ 10 times. 



(see Table II), therefore no special spectral feature 
seems to be visible. 

Concerning the electron spectrum produced by SN 
v e in a LAr TPC, in NH and sin 2 6 13 < 10~ 5 , 
we observe a bimodal distribution produced by a 
superposition of and F° x for E < -E sp iit and by 
mostly F® at high energies (see Table I) , producing 
a broad "shoulder" . Instead, in NH and sin 2 #13 > 
10~ 3 the electron spectrum will be mostly produced 
by the F^ . This is qualitatively similar to what 
happens also in IH. 

The presence of bimodal distributions and broad 
"shoulder" features in the spectra produced by v e 
and T> e for the Ci primary flux, then hold the 
promise of being observable at the detectors, in 
spite of the energy smearing. 

• In C2 the average energies of D e and v x are similar, 
as a result only one peak would be observed; the 



shoulder may not be discernible. Similarly in C3 
and C4, since the number fluxes of electron and 
non-electron flavors are close to each other, direct 
observation of the spectral split is hard. 

In summary, the spectral splits arc directly identifiable 
only in the phase C, when the average energy and 
luminosity of non-electron fluxes are sufficiently large. 

B. Earth matter effect 

As discussed in Sec. IIVB1 the passage of neutrinos 
through the Earth before reaching the detector can give 
rise to Earth matter effect oscillations in the spectra. For 
definiteness, in the following we will assume the neutrino 
path length in the Earth to be L = 8000 km. 

From Eqs. (|25p and (|26p. these oscillations have fre- 
quencies that are functions only of the solar neutrino 
mixing parameters, and hence are well known. It is then 
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FIG. 4: The same as Fig. [3] but for the benchmark flux model Ci. 



simply a question of identifying the oscillations in the 
final spectra, if they exist. As can be seen in Figs. [3] 
and 21 scintillation detectors and LAr TPCs can allow 
us to detect these spectral modulations. The following 
statements can be made from Table I and II and from 
Figures 3 and 4: 

• For the phase A, oscillations are expected in the 
v e spectrum for NH with sin 2 13 < 10~ 5 . These 
oscillations should be clearly detectable in LAr 
TPC. For NH with sin 2 13 > 1CT 3 , as well as for 
IH, there are no expected Earth effects in the v e 
spectrum. 

• For the phase A, oscillations are expected in the 
v e spectrum for IH with sin 2 6*13 > 10~ 3 and in 
NH. However in this case, the small flux differences 
between original antincutrino species leads to the 
presence of Earth induced oscillations at high ener- 
gies being barely visible, as can be seen from Fig. 3. 
The mixing scenarios IH with sin 2 6*13 < 10 -5 will 
not produce any Earth effects. 

• For Ci, Earth effects in v e spectrum are expected 
in all mixing scenarios (Sec Table I). In particular, 



for NH with sin 2 9 13 < 10~ 5 , the Earth effects are 
only at intermediate energies (10-25 MeV), while 
for the other mixing scenarios, they are prominent 
at high energies (E > 25 MeV). The sign of these 
effects is negative in NH and positive in IH. In 
particular, as one can see from Fig. 4, that the 
Earth matter modulations are clearly visible in the 
electron spectrum at high energies. This is due to 
the fact that in this case the spectral differences 
between F° and F° are relatively large at high 
energies. 

• For the case Ci, Earth effects in v e are prominent 
at high energies if the hierarchy is IH. For NH, 
the effects are expected at intermediate energies. 
However, due to the smaller spectral differences 
between F% and F% these features arc difficult to 
be observed, as can be seen from the Figure. 

• For C2, C3 and C4 - not shown in Figs. [3] and 
|4] - the Earth matter effects are in general less 
pronounced than in Ci. This is a result of some of 
the swaps being absent or only partially adiabatic. 
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the oscillations, one may rely on comparison between sig- 
nals at two detectors. It is quite possible that we will have 
more than one large water Cherenkov detectors within a 
couple of decades. Moreover, IceCube can be used as an 
accurate luminosity calibrator, and comparing the time 
dependence of luminosity in IceCube with that in another 
large water Cherenkov can help us ident ify E arth effects 
even without resolving the oscillations [1031 ] . Figure [5] 
shows some scenarios where the Earth effects may be 
observable without identification of oscillations, through 
the comparison between two detectors or the comparison 
between two different signals at the same detector. A 
combination of two 0.4 Mt-class water Cherenkov de- 
tectors, one shadowed by the Earth and the other un- 
shadowed [25| , has the capability to distinguish the mass 
hierarchies in the phase A for sin 2 #13 < 10 -5 , through 
an overall suppression of the event rate at the shadowed 
detector for NH (49|. In phase C, this combination can 
identify the hierarchies through the enhancement of high- 
energy events at the shadowed detector. Note that one 
of the two detectors could in fact be IceCube, as we do 
not really need the energy dependence. 



FIG. 5: Fhe left panel shows the relative asymmetry of 
events with and without Earth matter effects at two 0.4 Mt 
water Cherenkov detectors, for the benchmark flux model A. 
For Earth effects, the distance travelled through the Earth is 
taken to be 8000 km. In the right panel we show the same 
ratio for the benchmark flux model Ci. 

The above observations are consistent with the predic- 
tion that the Earth effects would appear only at those 
energies where the survival probabilities p or p in Tables 
HI and HIl are nonzero. 

The Earth induced modulations can be identified 
at a single scintillation or LAr TPC detector by the 
Fourier Transform technique [HI, IIH- Clear observa- 
tion of Earth effect oscillations in v e would then corre- 
spond to the possibilities (i) primary fluxes A, and NH 
with sin 2 #13 < 10~ 5 (ii) primary fluxes C. Observation 
of Earth effect oscillations in v e flux at the shadowed 
detector is more likely for primary fluxes C, but may be 
possible for A as well. The combination of the Earth 
effect signals in v e and v e may then enable us to shortlist 
the possible combinations of the flux and mixing scenar- 
ios. The Earth effects are visible in the intermediate 
or high energy range depending on the flux and mixing 
scenario. If one is able to divide a spectrum into energy 
regions where Earth effects are present and where they 
arc not, further distinction between flux and mixing 
scenarios is possible. However practically this task looks 
very hard. Also, note that non-detection of the Earth 
effects may be either due to their non-existence or simply 
due to their being too small to be detected. Therefore, 
in general the non-observation of these oscillations would 
not give us any concrete information. 

On the other hand, at a water Cherenkov detector 
whose energy resolution may be insufficient for discerning 



C. Suppression of v e in the neutronization burst 

The primary signal during the early neutronization 
burst (t < 20 ms) is pure v e . Since the model predictions 
for the ene rgy and luminosity of the burst are fairly 
robust [l04j . the observation of the burst signal gives 
direct information about the survival probability of v e . 
This probability is 0(6j 3 ) in NH with sin 2 13 > 10~ 3 , 
and sin 2 12 in all the other scenarios [l{|. Thus, the 
strong suppression of v e burst would be a smoking gun 
signal for the former scenario. Note that since no v are 
involved, self-induced oscillations are not developed and 
hence the collective effects do not give rise to any flavor 
transformations during the neutronization burst. 

However, in low-mass O-Ne-Mg supernovae (M ~ 
8 — 10 A/0), the MSW resonances may lie deep inside 
the collective regions during the neutronization burst, 
when the neutrino luminosity is even higher. In such 
a situation, neutrinos of all energies undergo the MSW 
resonances together, with the same adiabaticity |l05j . 
As long as this adiabaticity is nontrivial, one gets the 
"MSW-prepared spectral splits", two for normal hier- 
archy and one for inverted hierarchy [43|, 0, IH, [ill- 
The positions of the splits can be predicted from the 
primary spectra |44| . The splits imply v e suppression 
that is stepwise in energy. Such a signature may even be 
used to identify the O-Ne-Mg supernova, in addition to 
identifying the hierarchy. A LAr TPC with good timing 
and energy resolution would play a crucial role in this. 
However, we leave the study of this particular case for a 
future work. 
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D. Shock- wave effects 

Observables like the number of events, average en- 
ergy, or the width of the spectrum may display dips 
or peaks for short time intervals, while the shock wave 
is passing through the 7J-resonance. The positions of 
the dips or peaks in the number of events at different 
neutrino energies would also allow one to trace the shock 
propagation while the shock is in the mantle, around 
densities of p ~ 10 3 g/cc 0, |Tl| . The interplay between 
collective oscillations and shock- wave signatures has been 
recently studied in |6lj . This information, by itself 
or in combination with the corresponding gravitational 
wave signal, will yield valuable information about the SN 
explosion, in addition to confirming sin 2 #13 > 10 -3 and 
identifying NH (if the effects are seen in the v e spectrum) 
or IH (if the effects are seen in the v e spectrum) as the 
actual hierarchy. A caveat in this context is the role 
of matter turbulences in the post-shock regions and their 
impact on neutrino flavor conversions. Presumably, these 
would at least partially erase the signatures of the shock- 
waves in the neutrino signal [ill Ha - Hil ] . 

VII. SUMMARY AND CONCLUSIONS 

In this paper we have performed the first detailed study 
of the impact of collective and matter-induced flavor 
oscillations in the interpretation of the observable SN 
neutrino signal at a large water Cherenkov detector, a 
scintillation detector, and a Liquid Argon Time Projec- 
tion Chamber. We have analyzed the neutrino flavor 
evolution, including collective effects that give rise to 
spectral swaps and the MSW effects that further mix 
the neutrino flavors. In particular, we have taken into 
account the possible qualitative change in the fluxes and 
oscillation physics during the SN neutrino emission. We 
have calculated the neutrino flux at Earth, including 
possible Earth matter effects introduced if the neutrinos 
pass through the Earth before reaching the detector. We 
have also calculated the v e and v e spectra through the 
dominant channels at the above detectors, taking care of 
the detector efficiencies and energy resolutions. 

The collective effects, which are dominant for r & 500 
km, leave their imprints in the form of spectral swaps. 
The boundaries of spectral swaps are spectral splits, 
where the v e or v e spectra can have sharp jumps at 
critical energy values. Since the collective effects are 
essentially nonlinear, the number and position of these 
splits depends on the primary neutrino spectra. To 
explore this, we have fixed the average energies for 
v e and v e fluxes, assuming that their luminosities are 
almost equal (a result borne out by many SN neutrino 
simulations), and we have scanned the parameter space 
in the average energy and luminosity of the non-electron 
neutrinos. It turns out that, depending on the number 
and nature of spectral swaps, the parameter space may be 
divided into two main "phases": the phase A [L Va < L Vc ) 



that is typical of the fluxes during accretion, and the 
phase C (L„ x > L„ e ) that is typical of the fluxes during 
cooling. Phase A shows no spectral swap in NH while 
in IH shows a spectral swap at intermediate (10 MeV 
< E < 25 MeV) and high (E > 25 MeV) energies for 
both v and v. Phase C, on the other hand, shows spectral 
swap in both v as well as D at high energies even in NH. 
In IH it in general results in two spectral swaps, a e <H- y 
swap at intermediate energies and a e •<->■ x swap at high 
energies. 

We have also taken into account the MSW and Earth 
matter effects that determine the further neutrino flavor 
conversions after the collective effects are over. Armed 
with these results, we have looked for distinctive signa- 
tures of neutrino mixing pattern in the v e and v e spectra 
at the detectors, and we have examined the feasibility of 
their observation. In particular, we have considered (i) 
the sharp change in the spectrum that is the signature 
of a spectral split (ii) the Earth matter effects, (iii) the 
flavor conversion effects on the prompt v e neutronization 
burst, and (iv) the shock- wave effects in the neutrino as 
well as antineutrino spectra. 

The spectral split itself will be visible only if the cor- 
responding survival probability p or p changes suddenly 
at an energy > 10 MeV. While the splits in phase A only 
occur at low energies, the splits in phase C can occur 
at intermediate energies and can be observable. These 
splits will be a clear signature of the collective effects 
taking place deep inside the star, since no other known 
phenomenon can give rise to such a sharp change in the 
neutrino spectrum. Typically, the presence of spectral 
splits would produce bimodal observable energy spectral 
with two peaks corresponding to the ones of the initial 
v e and v x spectra. Our results for the other observables 
are summarized in Tables IIIII and IIVI 

From the tables, it can be observed that 

• Whether the Earth effects are visible at intermedi- 
ate energies, high energies, or not at all, depends on 
the neutrino mass hierarchy and the range of #13 . A 
clear identification of the Earth effects is therefore 
crucial in extracting neutrino mixing information 
from the observed spectra. (Note that since the 
threshold of the detectors will be ~ 5 MeV, we 
have ignored any effects at low energies.) If the 
Earth effects are present in the intermediate energy 
range and absent in the high energy range (or vice 
versa), it will be a signature of a spectral split. 
Thus, an indirect observation of a spectral split is 
also possible through Earth matter effects. 

• The vanishing of v & burst and the observation of 
shock wave effects is independent of the collective 
effects. These observables thus directy probe the 
neutrino mixing pattern: mass hierarchy and #13 
range. 

Concerning our study, we would to remind that while 
we have modeled the SN fluxes, and treated their oscil- 
lations in more detail than have been done in previous 
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TABLE III: Observable effects in the u e spectra, for fluxes of type A and C. Earth effects at low energies (E < 10 MeV) are 
almost impossible to detect and have not been considered. 





Phase A {L Ve > L Vx ) 
u e burst Earth effects Shock effects 


Phase C {L Ue > L Vx ) 
Vc burst Earth effects Shock effects 


NH sin 2 0i3 > 10" 3 
sin 2 0i3 < 10 -5 


Vanishes Absent Possible 
Present All E Absent 


Vanishes Only high E Possible 
Present Only intermediate E Absent 


ffl sin 2 0i3 > 10" 3 
sin 2 0i3 < 10 -5 


Present Absent Absent 
Present Absent Absent 


Present Only high E Absent 
Present Only high E Absent 



TABLE IV: Observable effects in the v e spectra, for fluxes of type A and C. The exceptions in region C are: with fluxes in 
regions C3 and C4 and IH. Earth effects are absent for sin 2 0i3 > 10 -3 , while they are present at intermediate as well as high 
energies for sin 2 0i3 < 10 -5 . Earth effects at low energies (E < 10 MeV) are almost impossible to detect and have not been 
considered. 





Phase A (L„ e > L Vx ) 
Earth effects Shock effects 


Phase C > L Vx ) 
Earth effects Shock effects 


NH sin 2 13 > 10" 3 
sin 2 13 < 10" 5 


All E Absent 
All E Absent 


Only intermediate E Absent 
Only intermediate E Absent 


m sin 2 13 > 10" 3 
sin 2 0i3 < 10" 5 


Intermediate and high E Possible 
Absent Absent 


Intermediate and high E Possible 
Only high E Absent 



literature, the results shown here are for a simplified 
treatment of the SN neutrino problem which ignores pos- 
sible trajectory-dependent effects, inhomogeneities, etc. 
All these effects would potentially introduce additional 
layers of complications in the simulation of the super- 
nova neutrino signal and their impact will need careful 
dedicated studies. In this sense, we believe that the 
predictions of possible signatures of supernova neutrino 
oscillations in large underground detectors need further 
investigations. This task is particularly timely now when 
different classes of large underground detectors are cur- 
rently under study for low-energy neutrino astrophysics 
and long-baseline experiments. In particular, from the 
perspective of supernova neutrinos it would be extremely 
useful to have as many different detection techniques as 
possible. 

We would like to emphasize that the information 
obtained from the SN v e spectrum is as crucial as 
that obtained from the v e spectrum, and in some cases 
even more useful. For example, the difference between 
the fluxes F® c and is always much more than the 
difference between the fluxes i 7 ^ and F^. Therefore, 
the features of flavor transformations like the spectral 
splits and Earth effects are likely to be more prominent 
in the v e spectrum, if present. Moreover, the shock-wave 
effects can be visible only in either v e or v e channel, 
depending on the mass hierarchy. The information from 
the v e spectrum observed at a LAr TPC would therefore 
be not just complementary to the one obtained from a 
water Cherenkov or scintillation detector, but it will also 
probe features of the SN neutrino signal that are not 
accessible to a v e detector. 



We wish to stress that the physics potential of large 
neutrino detectors proposed for low-energy neutrino as- 
trophysics is immense, and can be exploited further for 
studies of SN neutrinos. The accurate timing information 
obtained from water Cherenkov, scintillation and LAr 
TPC detectors should also allow us to track the SN v 
light curve quite faithfully. A detailed study of the time 
evolution of the neutrino signal could offer new insights. 
A LAr TPC would also be the most efficient detector 
for observing the v e neutronization burst. A comparison 
between the onset t iming of this signal with the onset 
time signals from D e |l06l - ll09l | detectable in water (or ice) 
Cherenkov and scintillation detectors, and gravitational 
wave signals may shed more light on the SN dynamics 
and neutrino emission. 

In conclusion, large future neutrino detectors would 
offer unprecedented opportunities to study supernova 
neutrinos and determine fundamental neutrino proper- 
ties through high-statistics studies of energy and time 
spectra. With the complementary physics potential of 
future water Cherenkov, scintillation and liquid Argon 
detectors, they promise to advance our understanding of 
the physics of neutrinos, and of their flavor conversions 
during a stellar collapse. 
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